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Abstract 

Lattice regularization of chiral fermions is an important development of the the- 
ory of elementary particles. Nontheless, brute force computer simulations are very 
expensive, if not prohibitive. In this letter I exploit the non-interacting character 
of the lattice theory in the flavor space and propose a multigrid approach for the 
simulation of the theory. Already a two-grid algorithm saves an order of magnitude 
of computer time for fermion propagator calculations. 
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1 Introduction 

After many years of research in lattice QCD, it was possible to formulate QCD with chiral 
fermions on the lattice jl], 0, ^ [|. 

The basic idea is an expanded flavor space which may be seen as an extra dimension 
with left and right handed fermions denned in the two opposite boundaries or walls. 
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Let N be the size of the extra dimension, Dy/ the Wilson-Dirac operator, and m the 
bare fermion mass. Then, the theory with Domain Wall fermions is defined by the action 

S DW := = Eli M(D 11 ~ l)A + P+ipi+i + P_Vi-i] (1) 

P+(i/; N+1 + mipi) = 0, P_Oo + mip N ) = 

where Ai is the five-dimensional fermion matrix of the regularized theory and D" = 
M — Dy/ with Me (0,2) being a mass parameter. 

I define also a theory with Truncated Overlap Fermions in complete analogy with the 
domain wall fermions by substituting 0: 

P+ -(£>" + 1)P+, P- - (D 11 + 1)P_ (2) 

Both theories can be compactified in the walls of the extra fifth dimension as low 
energy effective theories given by the action: 

S eff = XiDxi ^ 

xi = P+i'i + p~4>n, x\ = 4>%P- + ^nP+ 

where D is the chiral Dirac operator satisfying the Ginsparg- Wilson relation 0: 

75/r 1 + zrS = «7sP, (4) 

where a is the lattice spacing and R is a local operator trivial in the Dirac space. 

I defined Truncated Overlap Fermions such that in the large N limit one obtains 
Overlap Fermions || with the Dirac operator given by 0]: 

1 + m l-m 

Dov = — — 75Sgn(#) (5) 

where H = 75-D". 

Untill now, the computations with chiral fermions with the standard algorithms have 
been very expensive. The extra fermion flavors introduce a large overhead. One multi- 
plication with the fermion matrix costs 0(n) Dvy-rriultiplications with n ~ N for domain 
wall fermions and much larger for the overlap operator [Rl BL |T0f. 
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In this letter I propose a multigrid algorithm along the fifth dimension which makes 
these simulations much faster. The key observation is the lack of gauge connections along 
this dimension. It is well-known that the overhead of such algorithms scales like iVlogiV. 

Here it is the multigrid algorithm: ALGORITHMl (Generic) for solving the system 
D ov x = b: 

Given N, xq, r (= b), tol, toh, set tola = 1 and iterate: 
for i = 1, . . . 
tol = tol Q toli 

Solve Dy = r;_i within tolo (6) 
Xi = Xi-i + y 
n = b- DovXi 
ifWiWi < tol, end for 

where by o is denoted a vector with zero entries and toll, tol are tolerances. tol\ is typically 
orders of magnitude larger than tol such that the work per Dqv inversion is minimized. 

Remark 1. The straightforward application of the ALGORITHM! gives a two-grid 
algorithm. By calling it again in solving the smaller system and iterating, one gets a full 
multigrid algorithm. 

Remark 2. The corresponding Hybrid Monte Carlo (HMC) algorithm can be obtained 
by working with an approximate Hamiltonian in the coarse lattice and by a global cor- 
rection on the fine lattice. 

In Fig. 1 we compare the norm of the residual = b—Dov%i of the Conjugate Residual 
(CR) algorithm (which is optimal since Dov is normal |TTJ) and ALGORITHMl. I gain 
about an order of magnitude (in average) on 30 4 4 -configurations at /3 = 6.0 and m = 0.1. 

For the coarse lattice we used N = 6 with the Truncated Overlap Fermions and the 
Lanczos method to compute Dqv @- 

Let P be the following orthogonal transformation: 

Xl = P+ipi + P-ipi-x (7) 
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I computed the inverse of D in the ALGO RIT H Ml using 

IT 1 = (P T M- 1 M 1 P) 1>1 (8) 

where TWi is the same as M. but with m — 1 [Bj and the subscript 1, 1 stands for the 
(1,1) block of an NxN partitioned matrix along the fifth dimension. 

Recently, the possibility of a Multigrid algorithm along the all dimensions is raised 
I2fl . In this case a gauge fixing is needed. 

I would like to thank Philippe de Forcrand for suggestions on how to improve the 
ALGORITHM1 and Herbert Neuberger for interesting discussions after my talk at Lat- 
tice99 conference. 

The author thanks PSI where this work was done and SCSC Manno for the allocation 
of computer time on the NEC SX4. 

References 

[1] D.B. Kaplan, Phys. Lett. B 228 (1992) 342. 

[2] Y. Shamir, Nucl. Phys. B 406 (1993) 90; V. Furman and Y. Shamir, Nucl. Phys. B 
(1995) 54. 

[3] R. Narayanan, H. Neuberger, Phys. Lett. B 302 (1993) 62, Nucl. Phys. B 443 (1995) 
305. 

[4] P. Hasenfratz, V. Laliena and F. Niedermayer, Phys. Lett. B 427 (1998) 125. 

[5] A. Borigi, Truncated Overlap Fermions, Talk at The XVII International Symposium 
in Lattice Field Theory, Pisa, June 29 - July 3, 1999. 

[6] P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25 (1982) 2649. 

[7] H. Neuberger, Phys. Lett. B 417 (1998) 141, Phys. Rev. D 57 (1998) 5417. 

[8] A. Borigi, Phys. Lett. B 453 (1999) 46. 



[9] H. Neuberger, Phys. Rev. Lett. 81 (1998) 4060. 

[10] R. G. Edwards, U. M. Heller and R. Narayanan, FSU-SCRI-98-71, and 
[hep-lat/9807017 . 



[11] A. Borigi, Krylov Subspace Methods m Lattice QCD, PhD Thesis, CSCS TR-96-27, 
ETH Zurich 1996. 



[12] , L. Giusti, Ch. Hoelbling and C. Rebbi, |hep-lat/9906004 



5 




=E=10 



3 4 5 

D w multiplications 



x 10 



Figure 1: Norm of the residual error vs. the number of D w multiplications on 30 con- 
figurations. Circles stand for the straightforward inversion with CR and stars for the 
ALGORITHM!. 
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